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' Abstract. In this paper, we consider Beta(2 — a, a) (with 1 < a < 2) and related A- 

coalescents. If denotes the length of an external branch of the n-coalescent, we prove 
the convergence of n a ~ 1 T^ > when n tends to oo, and give the limit. To this aim, we 

. give asymptotics for the number of collisions which occur in the n-coalescent until the 

end of the chosen external branch, and for the block counting process associated with the 
n-coalescent. 

o\ : 
t— i . 

p/j ■ 1. Introduction 

1.1. Motivation and main results. In modern genetics, it is possible to sequence whole 
■ genomes of individuals. In order to put this information to maximal use, it is important 

to have well- fitting models for the gene genealogy of a sample of individuals. The standard 
model for a gene genealogy of a sample of n individuals is Kingman's n-coalescent (see [26], 
}27|). Kingman's n-coalescent is a continuous-time Markov process with state space the 
| set of partitions of {1, . . . , n}. The process starts in the trivial partition ({1}, . . . , {n}) and 

£> ' transitions are only possible as mergers of exactly two blocks of the current state. Each such 

binary merger occurs with rate 1. These mergers are also called collisions. 

For many populations, Kingman's n-coalescent describes the genealogy quite well. King- 
(T) ■ man showed in |27] that the ancestral trees of a sample of size n in populations with size 

, N evolving by a Wright-Fisher model will converge weakly to Kingman's n-coalescent for 

N — > oo (after a suitable time-change). This result is relatively robust if population evo- 
lution deviates from the Wright-Fisher model (see [27] or |28j). However, there is evidence 
that there are populations where the gene genealogy of a sample is not described well by 
Kingman's n-coalescent. Examples of such populations can be found in maritime species, 
^ , where one individual can have a huge number of offspring with non-negligible probability 



^ ! (see PQ, p2] [22], [20] and [IS]). 

A whole class of potential models for the gene genealogy of a sample was introduced inde- 
pendently by Pitman and Sagitov (see [30] and [31]): The class of n-coalescents with multiple 
collisions. A n-coalescent with multiple collisions is a continuous-time Markov process with 
state space *p( n \ where all possible transitions are done by merging two or more blocks 
of the current state into one new block. Every n-coalescent 11^ is exchangeable, meaning 
ToIl( Tl ) = n( n ) for every permutation r of {1, . . . , n}. The transition rate of a merger/collision 
of k of b present blocks is given by 

(1) \k= r x k (l-x) b ~ k x~ 2 A(dx) 
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for a finite measure A on [0, 1] (this definition is due to Pitman |30j). Since the process 
is characterized by the measure A, it is also called a A- n-coalescent. Note that Kingman's 
n-coalescent is a A-n-coalescent with A being the Dirac measure So in 0. 
It is possible to define a continuous-time process II with state space V, the partitions of the 
natural numbers N, whose restriction on {1, . . . , n} is a A-n-coalescent for all n £ N. Such a 
process is called a A-coalescent. 

An important subclass of A-n-coalescents are Beta n-coalescents characterized by A being a 
Beta distribution, especially for the choice of parameters 2 — a and a for a 6 (0,2). The 
class of Beta(2 — a, a)-n-coalescents appears as ancestral trees in various settings. They 
appear in the context of supercritical Galton- Watson processes (see [32] ) , of continuous-state 
branching processes (see [9]) and of continuous random trees (see [3]). They also seem to be 
a class where suitable models for the ancestral tree can be found for samples from species 
who do not fit well with the Kingman-setting (see [8]). Note that for a — > 2, the rates of the 
Beta{2 — a, a)-n-coalescent converge to the rates of Kingman's n-coalescent. In this sense, 
Kingman's n-coalescent can be seen as a border case of this class of Beta n-coalescents. 

For q = 1, Beta{2 — a, a) is the uniform distribution on [0,1]. The corresponding n- 
coalescent is the Bolthausen-Sznitman n-coalescent. It appears in the field of spin glasses 
(see [H], p2]) and is also connected to random recursive trees (see [2T]). 

Let us denote by IlM = (n| n) )j>o a n-coalescent. In this paper, we are interested in three 
functionals of n-coalescents 

• the length of a randomly chosen external branch ; 

• the number of collisions which occur in the n-coalescent until the end of a ran- 
domly chosen external branch ; 

• the block counting process 

R (n) = ( R (n) 

)t>o- = |n^ | is the number of blocks of 

Note that T^> can also be characterized as the waiting time for the first collision of a randomly 
chosen individual and o-W as the number of collisions we have to wait to see the randomly 
chosen individual merge. For i G {1, . . . , n} define 

lf> :=inf{t|{0^nW} 

as the length of the ith external branch and 

o-W :=inf{fc|{i}^7r fc } 

as the number of collisions until the end of the ith external branch, where is the state 
of the n-coalescent after k jumps. Due to the exchangeability of the n-coalescent, we have 

T (n) £ T (n) 

a( n ) = <j( n ) . Since we are only interested in distributional results, for the 

remainder of the article we will identify with and with a^. 

If the n-coalescent is used as a model for an ancestral tree of a sample of individuals/genes, 
the functionals and can be interpreted biologically. The length of an external branch 
measures the uniqueness of the individual linked to that branch compared to the sample, since 
it gives the time this individual has to evolve by mutations that do not affect the rest of the 
sample (see the introduction of [13] for more information). It was first introduced by Fu and 
Li in [18], where they compare mutations on external and internal branches of Kingman's 
n-coalescent in order to test for the neutrality of mutations. 
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The functional was first introduced in [T3], though n — a^ was also analyzed in [10] as 
the level of coalescence of the chosen individual with the rest of the sample. In both articles, 
the functionals were defined for Kingman's n-coalescent. 

For the biological interpretation of a^ n \ we see the n-coalescent as an ancestral tree of 
a sample of size n. Each collision in the n-coalescent then resembles the emergence of an 
ancestor of the sample, — 1 is the number of ancestors of the sample which emerge before 
the most recent ancestor of the randomly chosen individual/gene emerges (recall that time 
runs backwards in the n-coalescent). In this line of thought, gives the temporal position 
of the first ancestor of the chosen individual/gene among all ancestors of the sample, which 
are the collisions in the n-coalescent. Thus, gives the relative temporal position of 
the first ancestor of the chosen individual/gene among all ancestors of the sample (until the 

(n) 

most recent common ancestor). In this sense, we interpret as a measure of how ancient 
the chosen individual/gene is compared to the rest of the sample. 

In this article, we focus on the asymptotics of T^ n \ and for n — > oo. The 
asymptotics of these functionals are already known for some A-n-coalescents. For JW, we 
have 

• A = Sq (Kingman's coalescent): nT^ -4 T, where T has density t i— > M+t) 5 ( see 

m, eu), 

• A = Beta(l, 1) (Bolthausen-Sznitman coalescent): log(n)T^ n ) — > Exp(l) (see |17j). 

• A with a_ x = $ Q l x- l k{dx) < oo: 4 Exp(^{) (see [29], see also [15] ) 
for n — > oo. For a^ n \ we have 

• A = S : o-W /n A Beta{\, 2) (see [H]), 

• A = Beta{l, 1): 4 Beta{l, 1) (see [UJ), 

• A with n_ 2 = f^ x - 2 A{dx) < oo: trW A Geo(^) (see [19]) 

for n — > oo, where Geo(p) is the geometric distribution on N with parameter p. 

There are some known results for the asymptotics of the block counting process Define 
R = (Rt)t>o by setting Rt := |IIf |, the number of blocks of a A-coalescent II = (II^oo. Note 
that R = \\m. n ^ 00 R^ n \ if the block counting process R^ is defined for the restriction of IT 
to {1, . . . , n} for each n € N. For the small-time behaviour of the block counting process R, 
the following results are known: 

• If A(dx) = f(x)dx, where f(x) ~ Ax 1 ' 01 for 1 < a < 2 and ~ means that the ratio 
of the two sides tends to one as x — > 0+. Then 

(2) lim W-VRt = (— r) 1/(a_1) o..a. 

(see [3] and 0]) 

• If Jj e x j x~ 2 A(dx) = e~ a L e with e € (0, 1) and L e slowly varying as e — > 0. Denote by 
g t := L~ l e a ~ x , then 

\™ + eR gt = (r(2 - a)) 1/(1 " a) , in Probability. 

(see ®). 

Notice that both cases are complementary. The block counting process of Kingman 
coalescent gets similar limit value almost surely by taking a = 2, (see [3], page 216). 
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• For any A coalescent that comes down from infinity(that means for any t > 0, R4 < oo 
almost surely), there exists a deterministic positive function vt on (0, oo), such that: 

hm — = 1, a.s. 

t->-o+ v t 

For explicit form of v and finer results, we refer to [2]. This result is much more 
general than the two former ones. 

We will analyze the asymptotics for T", and for A-n-coalescents with A fulfilling 

pit) = C t~ a + 0(t~ a+ (), t -> 

for some Cq > 0, a £ (1, 2) and ( > 1 — 1/a, where pit) = L x~ 2 A(dx). Note that this class 
of n-coalescents includes all Beta(a, 6)-n-coalescents with parameters a € (0, 1) and b > 0. 
In this class of n-coalescents, we have the following asymptotics for T^ n \ and R^: 

o-(n) d 

• -r — tt — >■ cr, 

n(a—l) ' 

• for any t > 0, e > 0, P( sup In^i^Lc - (1 + C r(2 - a)t)" ll{a ~ l) \ > e) -> 0, 

0<t<i 

for n — >• oo , where a = Beta(l, a). 

Remark that if A(dx) = f(x)dx, where f(x) ~ Ax 1-01 for 1 < a < 2, we get that for fixed 
t > 0, 

iim i(i + ^^t)V(-^; : L = i 

in probability, which is to be compared to the following a.s. convergence implied by ([2]): 

hm 1( ^(2-^) ^ = L 

n->oo n a 

The different scaling comes from the fact that we do not consider the block counting 
process R = limj^oo R^ of a A-coalescent II, but instead the time-changed, scaled limit 
linin^oo To our best knowledge, there does not seem to be a direct link between 

these two results. 

Note that if we see the Bolthausen-Sznitman n-coalescent as a Beta{\, l)-coalescent and 
Kingman's n-coalescent as the borderline case of a Beta distribution with parameter a — > 2, 
the convergence results for shows a nice continuity in the parameters of the limit distri- 
butions in the range of Beta(2 — a, a)-n-coalescents with a S [1,2]. Our convergence result 
itself is even somewhat true in the border cases 1 and 2 (if one wages (a — 1) _1 — > oo for 
a — > 1 against log(n) — > oo for n — > oo for the Bolthausen-Sznitman n-coalescents). Also 
note that T obtained as the limit variable of nT^ in Kingman's case has the same law as 
2((l-cr)~ 1 -l), which gives again a nice continuity in results. The continuity also appears for 
the block counting process: replacing a by 2 in the formula for the Beta(2 — a, a)-n-coalescent 
gives the formula for Kingman's, this result will be proved in the sequel. 

Finally we can remark that together with the known asymptotics 4 a-1 for this 

class of n-coalescents (see [15], [20] and [23]), we have (jW/t^ -4 a. 

To prove these results, we will exploit some techniques from [15] . In [15] , they were used 
to analyze the asymptotics of a part of the length of a n-coalescent and the number of 
collisions in a n-coalescent for the same class of n-coalescents as analyzed in the present 
paper. For the convergence result for T^ n \ we present two proofs. One mimics the approach 
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in p3] and [23], using the convergence result for and = Yli=i Ti, where Tj is the 
waiting time between the (i — l)th and ith collision/jump of the n-coalescent. The other 
proof is based on the representation of as the first jump time of a Cox process driven 
by a random rate process which depends only on the block counting process associated with 
the remaining individuals labelled {2, 3, ... , n}. We use a recursive construction suitable for 
any A- n-coalescent: This construction consists in adding individual i to a coalescent process 
constructed by individuals from 1 to n except i such that consistence relationship is fulfilled. 

1.2. Organization of the paper. In section 2, we recall some known technical results which 
can all be found in [15]. In section 3, we obtain the asymptotic result about and also 
about the ratio between and r^ n \ Section 4 studies the small time behavior of the block 
counting process R^ n \ Depending on the property of R^ n \ our first method taking as 
the first jump time of a Cox process gives the asymptotic behavior of in section 5. In 
section 6, another method is provided by taking into account the fact that is the sum 
of a( n > initial waiting times for the coalescent process to jump from one state to the 
following. 

2. Preliminaries 



In this Section, we recall some results from |15j . 

Consider a n-coalescent with multiple collisions characterized by a finite measure A on 
[0,1]. Let v(dx) = x~ 2 A(dx) and p(t) = u[t,l]. When the process has k blocks, the next 
coalescence event comes at rate given by 



(3) 0k = £(/ + i) A M+i = J {o i f 1 " (1 " *)* ~ Ml " x) k ~ l 



\ A(dx) 
J x 2 



For n > 1, x G (0,1), let B nyX be a binomial r.v. with parameter (n,x). Recall that for 
1 < k < n, we have 

(4) mn„ > k) = [ t*-\i - *r* *■ 

Use the first equality in ([3]) and (jH) to get 

9n = j^{^jX k {l-x) n - k u{dx) 

= [ F(B njX >2)v(dx) 
Jo 

= n(n - 1) f (1 - t) n ~ 2 tp{t) dt. 
Jo 

All along this paper, the following hypothesis will be assumed 

(5) p{t) = C t~ a + 0(t- a+ <) 

for some Cq > 0, a € (1, 2) and ( > 1 — 1/a. Under this hypothesis, Lemma 2.2 of [15] gives 
us that, for n > 2, 

(6) g n = C r(2 - a)n a + 0(n a - min ^). 
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Recall that we call the number of coalescence events until reaching the common an- 

(n) 

cestor of the initial population (of size n). For k > 0, denote by Y" fc the number of blocks 
remaining after k jumps. Notice that is a decreasing Markov chain with Y Q = n and 
Yu = 1 for k > 7~( n ). Let -Xl = — Yu ^ e ^ ne num ber of blocks we lose during the 



kth coalescence event. We write X, 



(n) 




0. 



(n) 

The Markov property makes that the law of the first jump X\ will be of much interest. 

(n) 

We will look at some properties of X^ . Notice that 

(7) nx{ n) = k) = — F(B n , x = k + l)u{dx) 

9n Jo 

and thus 

(8) P(X (n) >k) = Io l nB n , x >k + l)u(dx) = (n-2)! - g^gp(g dt 

1 " 9n k\{n-k-l)\ £(l-ty-Hp{t)dt 

Under the same assumptions on p(t), setting eq > and 



(9) 



-c 


if 


c< 


a — 


1. 


1— a+eo 


if 


c = 


a — 


1, 


1-Q 


if 


c> 


a — 


1. 



Lemma 2.3 of [15] tells us there exists a constant Cfjo] s.t. for all n > 2, we have 



(10) 



E[X{ 



(«)l 



a — 1 



< 



Moreover, from Lemma 2.4 of [15], there exists a constant Cjjj] s.t. for all n > 2, we have 



E 



(») 



We consider n the Laplace transform of x[ n ^: for u > 0, (f> n (u) = E[e~ uX ^ ']. Assume 
that hypothesis © holds true. Let eo > 0. Recall (p n given by ([9]). Then we have (see [15], 
Lemma 2.5) , for n > 2, 



(12) 



4>n{u) 



+ 



7/' 



+ R(n,u), 



a — 1 a — 1 

where R(n,u) = (uip n + tt 2 ) h(n,u) with sup u6 [ ^] n>2 |/l(n, u)| < oo for all if > 0. 

Moreover, if we assume that £ > 1 — 1/a and set rj > — , then (from [15], Lemma 3.2) there 
exist £i > and Cjjj^iT) a finite constant such that for all n > 1 and u € [0, if], a.s. with 
a n = n~ v , 



Tn 



(13) 



^li^LI,^) <C^K)n-z\ 



We will also use the following result : Let V = (Vt,t > 0) be a a-stable Levy process with 
no positive jumps (see chap. VII in [5]) with Laplace exponent ip{u) = u a /(a — 1): for all 
u > 0, E[e~ uVt ] = e * u< V(<*-i). We assume that p{t) = C t~ a + 0{t~ a+ ^) for some C > and 
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£ > 1 — 1/a. Recall that 

T (n) is 

the number of coalescing events in the n-coalescent until 

reaching its absorbing state. Let 

|ntjAr(") 

V (n) =n -l/a J- (X^-^j) 
k=l 

for t G [0, a — 1), and 



r (n) 



-(») 



n ' | n 

fc=l 



Then, 

(14) (7 t (n) ,t6 [0,a-l])^(y t ,tG [0,a-l]) 

in the sense of convergence in law of the finite-dimensional marginals (see [15], Corollary 3.5, 
see also (201 125]). 

3. The number of collisions in an external branch 

Consider a n-coalescent with multiple collisions characterized by a finite measure A on 
[0,1]. Recall that u(dx) = x~ 2 A(dx) and p(t) = u[t, 1] which is assumed to satisfy ([5]). A 

n-coalescent takes its values in the set of partitions of {1, . . . , n}. For i > 0, let 7Tj = tt^ 
be the state of the process after the iih coalescence event. 

Pick at random an individual from the initial population and denote by the length of 
the external branch starting from it. Because of exchangeability, has the same law as 
the length of the external branch starting from the initial individual labelled by {1}. A 
quantity of interest will be a^ n \ the number of coalescence events we have to wait to see the 
randomly chosen external branch merging. Again because of exchangeability, cr^ has the 
same law as 

af ] =inf{i>0,{l} ^vr,J, 

the time we have to wait to see the external branch linked to individual 1 merging. We can 
write 

(n) 



(15) Zf > = £ 



-, 9 Y (n) 



where the e^'s are i.i.d. exponential random variables with mean 1. Note that the formula 
also holds true for and (just omit the subscripts). For the remainder of the chapter, 
we will identify with . 

In this section, we will determinate the asymptotic law of a^ n > for a class of coalescents 
containing the -Beta-coalescent with a € (1,2). 

Theorem 3.1. We assume that p(t) = C t" a + 0{t~ a+ () for some Cq > 0, a € (1,2) and 
( > 1 — 1/a. Then 

(!6) —. 7T ->■ o-, 

nya — 1) 

for n — > oo, where a = Beta(l, a). 
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Recall that in this class of n-coalescents, we also have r^/n -4 a — 1 (from [15] . see also 
[20] and [23] ) for n — > oo. Slutsky's theorem gives a convergence result for jW/tW, which 
measures how ancient the chosen individual is compared to the rest of the sample 

Corollary 3.2. We assume that p(t) = C t" a + 0{t- a+c >) for some C > 0, a £ (1,2) and 
£ > 1 — 1/a. TTien 



-(n) 



cr, 



/or n — >• oo, where a = Beta(l, a). 

proof of Theorem \3.1\ For convenience, we set 

(1) (*) = 0, if < a < 6, a G Z+, 6 € Z+, Z+ = {0, 1, 2, ...}. 

(2) log(0) = -oo. 

Notice that < r^ n \ Let y = (34> & > 0) denotes the filtration generated by Y^ n \ For 
any t > 0, we have 

P(o» > nt|3>) = P(cr ( " ) > nt,r (n) > nt|y) 

LntjAr( ?l ) 
i=l 



n 



NJ-W y (n) _ {X H + 1} 

11 v (n) 
i=l 



Notice that < 1 if i ^ r^, and = i jf j = r (»). 

We can hence write 



X— 1 i— 1 



LnijArW / (n) > 

log (P(o» > nt[y)J = ^ log 1 



x i-l 



and proceed to a power series expansion : 



log ( P(o» > n*|y) J = 4? + /, (2) 



with 
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, L«tjArM („) („) 

tf- E ocgci-^j + ^b+l). 

i=l r i-l r i-l 

where jj^ can be — oo if Xi ^ l = 1. Let us look further at I^j) . The idea is to replace 

(n) 

X by the limit of its expectation. 



j-(i) _ 7(1) , 7(2) 

1 nt ~ J nt J nt i 



with 

■£> = - E (^r+i)(^i 

8=1 

and 



a — 1 

We will use three lemmas whose proofs are given in the rest of the Section. 
Lemma 13.31 with rj = 1, tells us that, when < t < a — 1 

a — 1 Jo V a — 1 J \ a — 1 

Lemma 13.41 gives, for < t < a — 1, 

L«*JAt(") , , 

(18) 4? = -nVo-V-V- £ fen) _ 1 ) (yW) A 0. 

t=i ^ ' 
Finally, Lemma 13.51 gives, for t < a — 1, 

(19) 4? ^> 0. 
Adding dTTJ), ((THD and (US), we get that for i < a - 1 

log (p( ff W > nt|y)) 4 a log (l - , 

and thus 

P(a (n) > nt|y) 
While we know that P(cr (n) > nt\y) < l,then 



a — 1 



P(a (n) > raf) = E[P(<jW > nt|3>)] -> (\ 



a — 1 



We thus obtain that, for x € (0, 1), 

P(cr (n) > n(a - l)z) -»• (1 - x) a . 

a (n) 

and then that — converges in distribution to a Beta(l,a) law. □ 

n(a — 1) 
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Lemma 3.3. We set v v (t) = Jq (l — -^jj ^ dx,r] G R. We assume that p(t) = Cot~ a + 
0(t~ a+ {) for some C > 0, a G (1, 2) and £ > 1 - 1/a. For any < t < a - 1 and 

(1) Let to G [0, a — 1) and 5 > 0. XTie following convergence in probability holds when 
n — >■ oo: 

LntjAr( n ) 

n (a-l)/2-« sup 1^-1 V- M^V^^^I^O. 

0<*<*o V 7 

(2) Let t G [0, a — 1). T/ie following convergence in distribution holds when n — > oo: 
'( E MS) - n 1 '^)) ^ r, dr{l- r -)~^V r . 



Proof. The case n = a — 1 is given by Theorem 5.1 in [15] • Following the same arguments, it 
is easy to get the general result. □ 

Lemma 3.4. For any t < a — 1, the following convergence in distribution holds : 

\nt\ At'"' 

i=l ^ ' 

where (Vt)t>o is an a-stable Levy process with no positive jumps andv a (t) = J Q * ^1 — -^jj dx. 

Proof. Let 5 £ (0, a — I), to = a — I — 5 and t G [0, to]. 
Let e G (0, 1 — ■^ =l ) and /? = 1 - ^-y - e > 0. We have 



[ntjAr'") 

n i-i/a ^ (4 n, --^r)M ( - n iV = ^ + ^, 



i=i 



with 



and 



i=l 



1=1 V 7 

We will show that B n t converges to in probability and that A n t weakly converges to 
{v a {t)) l,a Vi as n — y oo. 

Convergence of A nt . Let z\ n ^ = n \ Y^{\ ^^ Y ^) >n py We have that sup nj j>! z[ n ^ < /3 _1 
a.s.. 

By using (JT2J) , it is enough to prove that 

E[exp(-nA ni )] ► e v«®* a /(a-l) 

n— ^oo 

for any n positive, where e v <x(t) ua /(a- 1 ) j s the Laplace transform of (ua(i)) 1 ' a 14- 
Taking uZ^ as Z- , we shall only consider the case u = 1. 
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Let us consider n^tf^ (4*)" = n^Efcf^ V/_1>^r We have > 

because the process {Y> n \i > 0) is decreasing, that 

i=i i=i 

P ({ y (Ln*jAr(«))-l < ^ 

(LntjAr("))-l 

i=i 

Use (|14p to get that the right-hand side of the last inequality converges to as n goes to 
infinity. Using also Lemma 13.31 with w = a, we have that 

|ntjAi-( n ) 

na -i ^ (yH)- a % a(f) , 
i=i 

as n — > oo. We can thus deduce that 

|nt|Ar( n ) 

(20) x; (^ B) ) B A« a (t), 

i=i 

as n — > oo. 

For a > 0, we set 



<2 = exp 



The process (M^l,k > 1) is a bounded martingale w.r.t. the filtration y. Notice that 



E[M^] = 1. As X? = and zf n) = for t > r( n ), we also have 



(fcAr< n ) , 
E (- 



fl- 1 /^^"' - log (n-^aZf } ) 



Using R(n,u) defined in (fl~2j) . we get that : 



M (a ] ,. 

n, [ntj 



exp 



i=l i=l i=l 



exp (-a^ ni ) exp ( -n" 1 £ V - \ [ - ]T R{Y^% n -^aZ^) 



i=l i=l 



Let 

^ a - 1 a - 1 



i=l 
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and write 

E [exp (—A„,t)] = A x + A 2 , 

with Ai=E [e~ Ant (l - e A ™)] and A 2 = E [e~ Ant e An ] . 

First of all, let us prove that A\ converges to when n tends to oo. Recall that the r.v 
are uniformly bounded by a.s.. Thanks to (|13|) . we have 

|ntjAr( n ) (2Z^\ a [nt\AT<- n '> 

1 E 



E 



M nlnt\ eX P I 



i=l 



a 



i=\ 



< M, 



where M is a finite constant which does not depend on n. By Cauchy-Schwarz' inequality, 
we get that 



(21) (A^ 2 < (E [e~ Ant |l — e An |]) 2 < E [e~ 



2A„t 



E 



;i- e ^) 2 



< ME 



'1-e 



The quantity A n is bounded and goes to in probability when n goes to infinity (see (|20l) 1. 
Therefore, the right-hand side of (|2ip converges to 0. This implies that lim n ^. 0O A\ = 0. 
Let us now consider the convergence of A%. Remark that 



A 2 = E 



\nt\ At' 71 ) 

fc=i 



a: 



Recall that E[M 



1. Using ([13]) . we get 

+ < A 2 < exp (q^/3- 1 : 



n Ci + 



! , Va(t) 



a — 1 



We get that lim„_ s . 0O A 2 = e Va ^^ a which achieves the proof. 

Convergence of i? n j. Here, we will use a similar approach as the one we used on the first 
half of pHU The process (Y^ n \i > 0) is decreasing. So if for some i < [nt\ , Y^\ < n(3, then 

(n) 

we have Y}'^ < np. Thus we get B nt = B nt l (n) 

1 (LntjArC"))-! Pl 



Moreover, 



(|ntjA-r< n ))-l 



(n) 



2=1 



a — 1 



)>*}, 



and then for any e' > 



(|ntjAr( n ))-l ' 

,(n) 



B n t\ > e') 



^({Y^^KnP) 

(|ntjA-r< n ))-l 



< P(n" 



E v. 



(n) 



5=1 



a — T 



>e). 



Use (|14p to get that the right-hand side of the last inequality converges to as n goes to 
infinity. 

□ 



(2) 

Now we deal with I nt . 
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Lemma 3.5. We assume that p(t) = Cot a + 0(t a+< ^) for some C > and ( > 1 — 1/a. 
Then, for any t < a — 1, we have 



(2), 

Tit I 



L"*JAt<") („) (n) 

E -o°g(i-^V- 



i=l 



V ( n ) 
z i-l 



Y; 



(n) 



i-1 



when n — > oo. 

Proof. Let < t < a — 1. First of all, remark that: 



0, 



(22) 
Indeed, 



Y, 



(») 



\nt\ p a — 1 — t 



11 



a - 1 



Y, 



(n) 



(LntJ)/(a-l) Et=i(4 n) 



n n n 

and we conclude using (fbl|) and the convergence of P(r( n ) > |_ n *J) to 1. Let us write 

irPOl 



nt 



with 



and 



,4 



n — \ ± nt 



/ (2) ll 



{l^j<(l-t/(a-l))n/2} 



- I^VjyW > ( i_ t/(a _i ))n/ 2}- 

The convergence (|22|) implies that A n tends to in probability. To prove the convergence of 
B n , let us first notice that for a € (0, 1), there exists a constant C{a) such that, if B n ^ x is a 
binomial r.v. with parameter (n,x), then 



(23) 



< 



E 



l{2<B n>x <(l-a)n}( lo g(l 



7? 



v{dx) < C(a). 



Indeed, there exists a constant C'(a) such that for u € (0, 1 — a), < — ln(l— it)— it < C'(a)u 2 . 
Hence, 



< 



E 



o 



1 {2<B Il , :E <(l-a)n}( 1 Og( 1 



11 



n 



u{dx) 



< C'(a) / E 

Jo 

< C'(a) / E 

J o 

<2C"(a) / E 
Jo 



B 



L {2<B n , x <(l-a)n}\—^-) 



v(dx) 



, B n ^ x 2 
n 



v{dx) 

Bn,x{Bn,x 1) 



11- 



v(dx) 



n- 



< 2C"(o)A([0,l]) =: C(o). 
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(2) 

Let us set a = (1 — t/(a — l))/2. Hence B n = \I^ t \lr (n) i . Notice that if n is large 

\ Y [nt\- n f 

enough such that an > 2, then if Y^™, > an we have > nt. Moreover, if Y^ n }, > an, for 

° — ' \nt\ — ' \nt\ — ' 

i < nt, we have > an > aY™ and = - Y^ < (1 - a)Y^\ < (1 - a/2)Y^\. 
Using ©, ([23]) and ©, we get that 

E[B n ] 

l nt \ y (n) , 1 x (n) i 

<^E[E[-(log(l - -L_^) + -^)l {1<Jf ,» <( , fl) ,. 1} l {y ( 1>m} l^]] 

i=l x i-l 

< C(q/2)nt ^ 

when n tends to oo. This achieves the proof of the Lemma. □ 

4. A RESULT ON SMALL-TIME BEHAVIOR OF THE BLOCK PROCESS 

We now turn to the study of the length of an external branch picked at random, denoted 
by . For any integer k between 1 and r^ n \ define as the time when the kth. jump 
is achieved. This variable can be expressed as a sum of k independent exponential random 
variables. More precisely, 

,(n) _ \ -« G.j 



4" = £ 



where the e$'s are independent standard exponential variables. Notice that = A^ n y We 

will first study asymptotics of A^\ For this, we use a two-step approximation method close 
to Section 4 of [15]. Define first 

aV = £ — > 



obtained replacing the e^'s by their mean, and 

1 

C r(2 - a) 



i(n) _ 1 ( V (n)\-a 

i=l 



obtained replacing by its equivalent in (0). 

Proposition 4.1. We assume that p(t) = Cot~ a +0(t~ a+( >) for some Co >0and(> l—l/a. 
Then, for any t < a — 1, we have 

when n — > oo. 

The proof is a straight consequence of Lemma 13.31 with rj = a and the following Lemmas 
^landOl 
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Lemma 4.2. Under the assumptions of Proposition \4^1\ we have 



n 



a-lf jW 



\nt\ 



0. 



when n — > oo. 
Proof. Use © to get 



.4 



(n) 



E 



Y; 



(n) 



i-1 



j=l 



- min(£,l) 



The result then follows from Lemma 13 .31 with r/ = a + min(£, 1). 
Lemma 4.3. Under the assumptions of Proposition f^TT] , we /iaue 

u>/ien n — >• cx). 



□ 



Proof. Recall that = (3^fcj k > 0) denotes the filtration generated by Y. Conditionally on 

e- — 1 

3^, the random variables are independent with zero mean. We deduce that 



9 v (n) 



E 



n 2a ~ 2 E 



sup 

t>0 



'[nt\AT^ 

E 

i=l 



< 4n 2a " 2 ^ 



i=i 



5y(") 



where we used Doob's inequality for the inequality. Thanks to ([6]) and Lemma 

/ \ 2 

,2a-2 ^LntjArW / 1 



= 2a, we get the 4n 2 " 2 X)i=i 



converges to in probability. 



with 



□ 



Heuristically, combining Theorem 13.11 and Proposition 14. 1\ we should get that n a ~ l T^ = 
n a-i^W^ conver g es j n i aw to (7 r(2- a ) (0- ~ (J ) 1 ~ a ~ !)• This line of proof will be followed 
in the last section. However, in the next section we will first present another way to prove 
this result with a method based on the consistency property of exchangeable coalescents. As 
a first step to this approach, we end this session with a result about small-time behavior of 
the block-counting process. 

Let denote the number of blocks of the n-coalescent at time t. The initial value 

Rq is n. We show that the limit law of the process R^ is deterministic under a certain 
time rescaling 

Theorem 4.4. We assume that p(t) = Cot~ a + 0(t~ a+< *) for some Co > and ( > 1 — 1/a. 
For any to > 0, e > , we have 

-(l + C r(2-a)t)- 1/(a " 1) | 



(24) 
when n — >• oo. 



( sup \n 1 R^l- 
o<t<t 



> e) -»• 0, 
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Proof. Let < r < a — 1, we have the following relation : 



A\ n) , L nr J ^ ■? 



Let i S [0, to], and define 
(25) r(t) = (q - 1)(1 - (1 + C r(2 - a)*)" 1 /^-!); 

on [0, to]. Notice that 

((!_ I^L)l-a_ 1)=t 



C r(2-a) vv a-1 

Then thanks to Proposition 14. 1| n a ~ 1 A^ r ^ t ^ converges in probability to t. 

Using the remark at the beginning of the proof in Lemma we get the convergence 



n-iijW = J^i (i _ I^L) = (i + CoT{2 - f 

^i„„^i i ^ a — 1 



\_nr(t)\ 

(n) 

when n — s> oo. Moreover, since R± is decreasing, then for any < 5 < 1, 
lim P(flW < i^L +1 < fiW ) = L 

The constant (5 being arbitrary, we thus obtain the convergence in probability of n~ 1 R^ 1 _ a 

to (l + C r(2-a)i)~ 1/(Q ~ 1) . 

We obtain (|24p using again the fact that R^ is a decreasing process. 

□ 

In fact, the asymptotic result concerning block counting process of Kingman's coalescent 
is also valid. The method is almost identical to that employed in the above Theorem. In the 
context of Kingman's coalescent, we use the same notations n( n ), A^ n , R^ n \ 

Theorem 4.5. In the setting of the Kingman's coalescent, for any to > 0, e > , we have 
(26) P( sup {n^R^U - (1 + V 2 )" 1 ! > e) -> 



0<t<t 



when n — >• oo. 



Remark that this Theorem shows a nice continuity from Beta-coalescent(the process that 
we consider is more general but contains Beta-coalescent) to Kingman's coalescent., setting 
a = 2 in Theorem 14^41 



Proof. Recall that is the time when ith jump is achieved. When Il( n ) has b individuals 
at some time t, then the process encounters the following coalescence at rate ( 2 ) where 
two randomly chosen individuals will be coalesced. Il( n ) remains 1 when all individuals are 
coalesced. 

For < t < 1, we have 



An) sr &k_ 

\nt\ Z_^i (k\ 

k=n-\nt\+l \2J 
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where ejS are i.i.d unit exponential variables. Notice that 

k=n-\nt\+l \2) L J 

as n tends to oo. There exist a constant K > 0, such that, 



fe=n-[n*J+l \2l 

So we deduce that 



nA \ n \ ,3 2(— - 1) = — := /(t) 

as n converges to oo. 

We denote by := f/(i + 2) the inverse function of /(t). 

Similarly, i?( n ) is decreasing, so 

F{B$> < 4% < R%> ) -> 1, 

v A lnf-Ht-t)Wi ~ tn A w-i (t+i )W] j 

as n tends to oo for any < 5 < t. 

So t"" 1 - L"/~ 1 ( t )J 4. o 
Furthermore, 

^.(n) y(») , ,_1/,m n 

= r Ln/-Mt)J = l n f (*)] ^ j _ = 1 

n n n 1 + 1/2 ' 

as n tends to oo. So — ^ >• 1+t / 2 • 

Using again the decreasing property of R\ , we finish the proof. □ 



5. The length of an external branch picked at random 

Dynamics of any exchangeable coalescent with multiple mergers are characterized by rates 
Afe k which suit a consistent relationship (this is Pitman's structure theorem, see |30j . Lemma 
18): 

(27) X^k = Ab+i^+i + Xb+i,k- 

This relationship comes from the fact that k given merging blocks among b can coalesce 
in two ways while revealing an extra block : either the coalescence event implies the extra 
block (and then k + 1 blocks will merge) or not. Thus we get a recursive construction of the 
ra-coalescent process U^ n \ 

Let us define IT n ' 2 ) as the coalescent process of individuals labelled from 2 to n. Now we 
consider the individual labelled by 1. The lineage of this individual can be 'connected' to 

n( n < 2 ) 

• either at any of its jump times, in which case block {1} participates to a multiple 
merger implying at least 3 blocks, and we call this collision "Type 1" (see Figure [1]), 

• or at any other time to one of the present blocks and then participates to a binary 
collision, and we call it "Type 2" (see Figure [2]). 
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1 . 



Figure 1. n = 5. Individual 1 is chosen. Type 1: individual 1 encounters a 
multiple collision. 



l , 



2«- 
3«- 



Figure 2. n = 5. Individual 1 is chosen. Type 2: individual 1 encounters a 
binary collision 



From now on, our analysis is conditional on U^ n ' 2 \ Between two jump times of H^ n ' 2 \ 
assuming that there are b blocks in Il( n ' 2 ), the extra block coalesces at rate b\ b +i.2- If the 
extra block remains unconnected just before a coalescence event implying k blocks among b, 
then it will participate to this event with probability 

fix k+1 (l-x) b - k v(dx)/g b _ i X b+lik 
flx k {l-x) b ~ k v(dx)/g b \k 



(28) 



This equality comes from (|27p . Let us see how to get the law of , the coalescence time of 
individual 1. We define by R^ n ' 2 ^ the block counting process of n( n ' 2 ). Notice that it has the 
same law as i?( n_1 ). We introduce 



T c (n) the first jump time of a Poisson process r]c directed by the measure u, 



,{n) 



K t A p (n,2) 



RV ' '+1,2 



dt: 



(n) 

the time of the first appearance of 'Head' in the following coin flip, independent 
of : at each jump time t of R,[ n ' 2 \ we toss a coin, and get 'Head' with probability 



fl("' 2 )+i,fl(^ 2 )-fl("' 2 )+i 

A J? („,2)^(„,2)_ fl (n,2) +1 



and 'Tail' with probability 



Then, conditionally on , and T c (n) A Tj n) have the same law 



i£' 2 )+i,«£' 2) -.R£ n ' 2 >+i 



,(n,2) (n,2)__(n,2) 



(see (EHD)- 



Remark 5.1. A more formal way to interpret 

T (n) 

is as follow. Let be Cox process directed 



by random measure + , where u^ 1 ' = ; 



{t is a jump time} A ( „, 2 ) D („,2) D (n,2) , , 
«t_ >" t _ -«i + 1 

5t is the Dirac measure in t (see |25(. p.226]). Then T^ n > has the same law as the first jump 
time of £( n ) 



Let us now give our main result 
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Theorem 5.2. The following convergence holds : 

d m 1 



"°" Ir< " ,Ar =cW2^)(( 1 -^ 1 "°- 1 »- 

for n — ^ oo . The density function of T is 

fT(t) = aCor{2 ; a \ l + C r(2-a)t)-^-\ t>0. 
a — 1 

In particular, in the Beta(2 — a, a) case, the density is 

Mt) = - — i— + )-^r- 1 , t>o. 

(a — ljr(a) aT(a) 

Proof. For the sake of simplicity, we will make the proof only in the Beta{2 — a, a) case. The 
proof can be extended to the more general case where ([5]) is satisfied with the details omitted 
here. In this special case, Co = (aT(a)r(2 — a)) -1 and dynamics are given by 

B(k — a,b — k + a) 



A 



b ' k B(a,2-a) 
where B{a, b) is a Beta function of parameters a and b. 

Define r^ n ' 2 ^ as the number of jumps of the process up to time n 1 a t. It is a 

straightforward consequence of Proposition 14.11 that 

(n,2) 

(29) - 1 ^ r(t), n-^oo 

n 

for t — y oo , where r(t) is defined in (|25p . 

For i > 0, in the process n^' 2 -*, we denote by Y^ n ' 2 ^ the number of blocks remaining after 
i jumps which equals 1 from the time all individuals are coalesced to 1, and Y^ n ' 2 ^ = n — 1. 
Let X^ n ' 2 ^ = — Y^ n ' 2 ^ be the number of blocks we lose during the ith. coalescent event. 

We write X<j n) = 0. Notice that (y("> 2 ), J£(™. 2 )) has the same law as ((F^" 1 ), X^" 1 ))) . 
Using the description given above, we have 
p / n a-i r (n) > ^ 



E[P(n a -'(T( n) ATj n) ) > t\U^)} 
E^n " 1 ^ > t|n( n ' 2 >)P(n a - 1 rj n) > t\n in '^)} 



rt rl , „ N _ , n . 2) . r ±_ A („. 2 ) 1 , „(n,2) 



E[exp(- / / ^^^^(l-x)^ 2 --" 1 !/^)^)!" 



(n) 

t R/'O _ ™ »( n ' 2 ) _L™_ r ' V ( n > 2 ) _ V-( n > 2 ) 



i=l 

We decompose the term in the expectation into two parts: the exponential on one side and 
the product on the other. 

Let us first look at the exponential term. Using Stirling's formula we get that, for < s < t, 

l-a R (n,2) B ( 2 ~ Q^ln'i 2 -" + a ~ 1 ) = 1-q (^in'i-^ / ^in' 1 -" \q-l f(f>M ) 

snl ~ a B(2-a,a) F(a) K n ' J{ snl ~ ah 
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where / = /(i){t>o} is a deterministic function which converges to as t converges to oo. 
The sequence (R^il a ,n > 2) is decreasing so, thanks to Theorem I4.4| we deduce that 

sup 0<s<f ( sn ^~ a ) a ~ l f {R^ n i- a ) converges in probability to as n tends to oo. Consequently, 
using again Theorem 14.41 we get that 

/"* , /_ 21 5(2 - a, i? (n \ 2) a + a - 1) P t 
(30) exp(- jf ^ Us) 4 (1 + _)-«, „ oo. 

Convergence of the product term is obtained by the same method as in proof of Theorem 
13.11 combined with the convergence in (|29|) . To avoid showing almost the same reasoning, we 
leave the details to readers. This way, we have 

" ( ™' 2) ,,(71,2) v (n,2) , , 



i=l -'i-l v 



The product of ([30]) and (|3ip then converges in probability to (1 + ^j^jy) ^- Since 

this product is bounded, we get that 

p/ n a-i T (n) >t \_>t 1 + n ^ oo. 

aT[a) 

We achieve the proof. □ 

As a consequence of Theorem 15.21 we can get an asymptotic result on the size of the 
population at the moment of collision of individual 1. 

Corollary 5.3. The following convergence holds : 

n " ly jS (1 + C T{2 - cOT)- 1 /^- 1 ) = 1 - a 

for n — > oo. Moreover, the density function of this limit is ax a_1 l|o<a;<i}- 

Proof. In terms of block counting process, we have Y^fy = Notice that R^( n ) = 

R^ n ) , ,,w„^. Using Theorem 15.21 we known that n a ~ l T^ converges in distribution to T. 
Hence, if to > 0, we deduce from Theorem 14.41 that 

w i-°.( n °.-ir(n)) d (T<rr(o ^tW/(«-i) 

This achieves the proof. □ 

6. AN ALTERNATIVE PROOF FOR THEOREM 15.21 

In this section, we present an alternative proof for Theorem 15.21 using the convergence 
results for from Theorem 13. 11 First, we need a stronger version of Proposition 14.11 which 
gives weak convergence in the path space. Recall that 



A [n) = V 

k h ^ 

where the e^'s are independent standard exponential variables. 
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Proposition 6.1. We assume that p(t) = CQt~ a +0(t~ a+( *) for some Cq > and C > 1— 1/a. 

Then, for any t < a — 1, we have 

(32) (n^A^ ) s<t A (—1 .((1 - - l)) s < t , 

L " SJ Cot(2 — a) a — 1 

m t/ie sense of convergence in the path space D[0,t] for n — > oo. 
Proof. Note that Theorem 14.11 states 

for < s < a — 1 and n->oo. So for every fixed s € [0,t], we have pointwise convergence in 
probability in (I32[> . This implies weak convergence of all finite dimensional distributions due 
to the subsequence criterion for weak convergence. In order to show weak convergence in the 
path space, we will show tightness for the distributions from (I32p . Since the limit process is 
continuous, it suffices to show that the condition (i) of Theorem 7.3] and condition (7.12) 
from Corollary 7.4] are fulfilled (see Corollary 13.4]). For the present processes, these 
conditions translate to showing that for every e > and ij > 0, 

(i) there exists a > s.t. P(n a ~ 1 A^ > a) < rj for n big enough and 

(ii) there exists a < 5 < 1 so that 

1 P f^a-lf An) An) 



y v Ln-min(ti+o,tn \n-t\\> — / — " 

for n big enough and any t\ G [0,t]. 
Condition (i) is trivially fulfilled, for condition (ii) we can use Theorem 14.11 to show that for 



n oo, 



Ln-min(ti+<5,f )J L n '*iJ 

where f(s) := Cor( ^_ tt) ((l - 7 ^ T ) 1 ~ a - 1). Note that P(/(min(t 1 + S,t)) - f(h) > e) < 
P{f(t) ~ fit ~ — e ) ^ {0' !}• Since / is continuous, you can now choose 5 small enough 
that f(t) — f(t — 5) < e and then n big enough to fulfill (ii). Thus, we have shown tightness 
of the distributions in (I32p which establishes the desired weak convergence □ 

Now we come to the alternative proof of Theorem 15.21 

Alternative proof of Theorem \5.SX Fix t G [0, a — 1). We have 



(a^/(n(a - 1)), (n"" 1 ^ ) s < t ) A (a, ( * _ ((1 - -^) l ~ a - 1)), 



<t 



a — l 

for n — > oo. Due to Skorohod-coupling, we can assume that this convergence also holds 
almost surely. Since s i-> ( Cor ^ 2 - a ) ((^ ~~ o^r) 1_a ~~ -0) ^ s continuous on [0, t], the almost sure 

convergence of (n a ~ 1 A^j) s <( in -D[0,i] is even almost sure uniform convergence on [0, t] (see 
p. 124]). For any series (rr n ) ne ^ on [0,t] with x n — > x, we thus have 

\nx n \ <7 r(2-a) U a-V ' 

almost surely for n — > oo. The only problem left is that a^ n \ a may take values in [0, a — 1) 
and not only in some subset [0, t]. To remedy this, note that if we restrict all random variables 
on {a < a — 1 — ^} for k E N, we have ^"'(uj/n < a — 1 — jr for n = re(u;) big enough 
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for almost all uj € {cr < a — 1 — V\. Thus, by using the Skorohod-coupling for the series 



(aW/lnia-l)),^- 1 ^)^^,) 



we have 



a-1 An) d + //-. U-q -, n 

V-M/nJ ^ C r(2-a) U a-l j J ' 

almost surely on {<r < a — 1 — i} for the coupled versions of these random variables (note 
that o» < r( n )). Since a is Beta-distributed, we have, for k — > oo, 

p({n a " 1 A( n) . . . . e •} n {a < a - l - £}) ~ FK-^J" 1 (n) , i G •). 

v L \n-a\ n > /raj ' L — j^. J 7 v p 1-0 ^ ^ /raj 7 

n ^ " n ^ C r(2-a) ((1 a-l j 1} - 

The only thing left to prove is that T has density fx- This is done by computing the 
distribution function 

■((i-r^T) 1_Q -i)<0 



This shows 



y C T{2-a) yy a-1 

P(<7 < 1 - (1 + C F(2 - £*)*)==*) 
l-(l+Cbr(2-a)t)S=T 

a(l — x) a dx 



= 1 - (1 + C r(2 - a)t)"=r. 
and finally by differentiating. □ 
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